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INTERNAL  TIDAL  WAVES  IN  THE  OCEAN 
By  B.  Haurwitz 

Abstract 

The  theory  of  internal  tidal  waves  is  briefly  reviewed  in 
f order  tc  demonstrate  the  importance  of  simultaneous  observations 

at  different  positions.  After  a description  of  the  expectancy 
test  it  is  shown  that  a reasonably  high  probability  exists  that 
the  internal  tidal  oscillations  found  at  Meteor  Station  335  in 
the  temperature  observations  are  real  and  not  purely  formal 
results  of  the  harmonic  analysis.  The  existence  of  internal 
tides  at  Meteor  Station  I4.38  and,  to  a lesser  degree,  at  an  anchor 
station  of  the  Altair  cannot  be  regarded  as  established  because 

i 

1 of  the  high  probability  that  amplitudes  of  the  magnitude  obtained 

i 

by  harmonic  analysis  would  occur  in  random  data.  The  study  indi- 
cates clearly  that  longer  series  of  observations  are  required  if' 
the  existence  of  internal  tidal  waves  is  be  established  securely. 
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1.  Introduction 

Internal  waves  :n  the  ocean  have  received  considerable  atten- 
tion in  recent  tears.  Apart  from  the  intrinsic  interest  whic: 
these  phenomena  ueserve  it  can  be  surmised  that  they  play  an 
appreciable  role  in  furthering  mixing  processes  in  those  regions 
of  the  ocean  where  they  occur.  This  has  become  especially  apparent 
after  Defant  (19U3)  has  made  it  very  probable  that  they  are,  at 
least  in  the  regions  of  the  continental  shell's,  suoject  to  break- 
ing, like  the  waves  at  „he  cc.can  surface  in  shoaling  water. 

For  the  practice  of  oceanographic  surveys  it  is  of  course 
necessary  uc  know  if  internal  waves  are  present  ir.  oraer  to  elimi- 
nate the  effects  of  such  waves  on  the  analysis  ana  interpretation 
of  the  data.  The  importance  of  this  correction  has  recently  been 
demonstrated  very  convincingly  by  Defant  (19p0). 

The  effect  of  waver,  at  the  ocean  surface  on  ship  motions  is 
well  known  and  cioes  not  need  ar.y  further  elaboration.  The  "uoad- 
water"  phenomenon  caused  by  internal  waves  generated  by  the  motion 
of  ships  has  also  beer,  fully  described  in  text  books,  for  instance 
by  Sverdrup  ar.d  Coll.  (19l;f).  in  a similar,  but  much  stronger 
way  sucmarlne  navigation  may  be  effected  by  internal  waves  at 
lower  depths.  And  it  appears  that  submarines  have  encountered 
such  internal  waves,  with  rather  unpleasant  results. 

As  far  as  the  excitation  of  internal  waves  is  concerned  it 
can  oe  surmised  that  they  are  in  part  due  to  atmospheric  causes, 
like  the  ocean  surface  waves.  Other  internal  waves  may  h*»  caused 
by  seismic  activity,  cut  in  view  of  the  relatively  infrequent 
occurr er.ee  of  earthquakes,  as  compared  to  the  apparently  quite 


frequent  occurrence  of  internal  waves,  it  is  necessary  to  attri- 
bute the  origin  of  at  least  many  of  the  observed  internal  waves 
to  other  causes. 

Since  one  of  the  most  important  forces  procucing  periodic 
changes  in  the  oceans  is  the  tidal  action  of  the  sun  and  the 
moon  particularly  attention  has  naturally  been  given  to  oscilla- 
ticns_ of  tiaal  serious  in  the  study  of  internal  waves. 

In  a previous  report  (Haurwitz,  1952)  various  papers  were 
listed  which  present  observations  of  internal  tices.  Further,  a 
series  of  temperature  soundings  made  by  the  Atlantis  from  16  to 
29  June  1936  at  latitude  3U°i«  and  6o°.v,  northwest  of  bermuaa  was 
analyzed  in  order  to  determine  the  existence  or1  nonexistence  of 
internal  tides.  A statistical  test  disclosed  that  no  internal 
waves  of  tidal  character  were  present  in  the  Atlantis  material, 
since  the  amplitudes  of  the  oscillations  determined  by  harmonic 
analysis  are  not  larger  than  can  be  expected  in  random  data. 

Tills  conclusion  suggested  further  chat  some  ox  one  ocher  reported 
occurrences  of  internal  tidal  waves  might  not  be  real,  but  simply 
due  to  the  formal  harmonic  analysis  of  the  observations. 

with  these  considerations  in  mind  some  further  ocservalion 
material  lias  since  been  subjected  to  the  same  statistical  tests. 
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flcally,  the  data  considered  are  those  obtained  by  the  German 
North  Atlantic  expedition  with  the  Meteor  in  193?  and  1936  v von 
Schubert,  19UU ) and  by  the  fnternat ional  Guifstream  Investigation 


the  Altair  1936 


;befant,  1940a)  which 


been  readily 


accessible  bo  the  author.  The  analysis  has  in  every  case  been 
confined  to  the  temperature,  as  in  tne  case  of  the  Atlantis  data, 
tor  these  three  series  of  observations  which  are  treatoa  here 
salinity  and  current  measurements  are  also  available.  but  these 
data  were  not  analyzed  since  the  necessity  of  statistical  tests 
can  be  demonstrated  quite  clearly  on  the  basis  of  the  temperature 
observations  alone.  Such  tests  require  a sufficient  amount  of 
data.  hence,  for  reasonably  reliable  determination  of  the  pre- 
sence of  internal  tides  and  their  nature  the  observations  should 
comprise  more  than  just  four  or  five  tidal  periods.  Apart  from 
this  suggestion  for  the  collection  of  further  data  on  internal 
tinal  waves  there  are  some  additional  data  which  are  indispensable 
if  a reasonably  coherent  picture  of  this  important  phenomenon  is 
to  be  obtained.  In  ox-der  to  explain  these  requirements  It  is 
necessary  to  consider  briefly  the  theory  of  internal  tidal  waves. 


x x , niom  ty  yjf  Liilcnml  ixues 

The  internal  wave  motions  which  are  possible  in  a stratified 
coan  have  been  studied  by  Kjeldstad  11933).  "be  has  also  given  a 
numerical  Integration  method  which  permits  the  computations  of 
the  possible  internal  waves  in  ar.  ocean  with  arbitrary  density 
d i ? t r 1 o'  t ion.  More  recently  Iroer.  ( 193^)  has  studied  internal 
waves  in  a fluid  with  a certain  specified  density  distribution; 
namely,  a rather  rapid  change  in  an  intermediate  liver  and  almost; 
constant  density  in  two  very  aeep  layers  above  and  below  the 
inhomogeneous  layer.  here  we  shall  confine  ourselves  to  the  dis- 
cussion of  a simple  case,  following  Fjeldstad's  work  without 
rev' owing  his  xa;er  in  detail. 


If  a rectangular  Cartesian  coordinate  system  is  chosen  with 
the  x anu  y axes  horizontal  ana  the  z axis  vertically  upward  the 
velocity  components  may  be  denoted  respectively  as  u,  v,  ana  w. 

Let  £ - ?_  (SO  1 r.  L he  the  Coriolis  parameter,  g the  acceleration 
of  gravity  and  ar.d  p the  density  and  the  pressure,  respectively. 
The  fluid  may  originally  be  in  a state  c!  rest  cn  which  the  inter- 
nal v;ave  motion  is  superimposed  as  a small  perturbation , Then  u, 
v,  w are  perturbation  quantities.  Further,  let  the  density  and 
pressure  in  the  undisturbed  state,  be  distinguished  by  the  subscript 
zero,  and  p0  while  the  total  density  and  pressure  in  the  dis- 

turbed case  are  ^ jQ  and  + jD . The  perturbation  quantities, 
u,  v,  w , p ,JD  are  assumed  to  be  sc  small  that  terms  of  higher 
order  can  be  neglected.  The  fluid  is  incompressible,  but  the  un- 
disturbed density  J3,  may  be  a function  of  z.  Then  the  three 
equations  of  motion,  the  equation  of  continuity,  ana  the  equation 
expressing  the  constancy  of  the  density  of  an  inuiviuuai  parcel, 
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The  dependence  of  these  variaole  on  x.  y ana  t may  oe  expressed 

by 
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This  expression  represents  a wave  progressing  in  the  positive  x 
direction  with  the  speed  c - 'y/Jz.  , lh>-  period  27 v <r  and  the  wave 
length  2 777-4,.  It  assumes  furthermore  an  exponential  or  periodic 
change  in  the  y direction  depending  or,  whether  the  factor  OC  is 
real  or  imaginary . 

If  all  the  other  unknown  functions  arc  expressed  by  w one 
finds  that 


where 


= JU 
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(2.3) 
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quations  (2.7)  ana  (?.&)  shew,  when  the  real  parts  ere  taken, 


that 
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Thus  the  hodograpn  cf  the  tidal  velocity  is  an  ellipse.  The  ratio 
cf  the  two  axes 

-»/>>  _ /-*-jF5t 


If  CV,  1 r complex  or  imaginary  the  foregoing  representation  of  the 
hodograph  would,  cf  course,  have  to  be  modified. 

The  ratio  of  the  two  semi -axes  can  be  obtained  from  observe- 
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enters  through  the  Coriolis  parameter  J.  the  wave  length  2 Tf/dt 
and  the  factor  oC  which  determines  the  dependence  on  y.  nence, 
even  if  the  perioc  is  given,  as  would  be  the  case  for  tidal  waves, 
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the  internal  wave,  if  it  it;  progressive,  and  would  thus  greatly 
increase  the  information  while  at  the  same  time  allowing  a more 
complete  check  on  the  theoretical  deductions.  Observations  by 
more  than  one  slip  have  been  made  by  rj elds tad  curing  a study  of 
internal  tidal  waves  in  a Norwegian  fjord  and  by  0.  Neumann  (1949) 
as  part  of  an  investigation  of  shorter-period  oscillations.  but 
from  the  open  ocean  such  data  have  net  ooen  available  at  the  time 
when  this  report  was  completed. 

The  differential  equation  (2.3)  can  be  simplified,  as  shown 
in  detail  by  Fj  elds  tad,  since  the  small  variation  of  the  undis- 
turbed density  may  be  neglected  in  the  first  term  on  the  left- 
hand  side,  so  that 
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*- 


(2.3a) 


In  order  to  make  the  croblem  determined  it  Is  necessary  to 
Introduce  boundary  conditions.  At  the  same  time  the  particular 
fluid  model  to  be  discussed  here  briefly  will  be  specified.  It 
may  consist  of  a homogeneous  layer,  1,  extending  from  z = 0 to 
z = hj  with  a rigid  lower  boundary,  and  of  an  upper  layer,  2, 
with  variable  density,  extending  from  z = hp  to  z = h = hp  * hp. 
the  free  upper  oounaary.  At  the  internal  boundary  the  density 


is  assumed  to  ce  continuous, 
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At  the  lower  ri^iv.  boundary  the  vortical  velocity  must  vanish, 
w = « ‘ a c 7,-0 

At  the  internal  boundary  the  vertical  velocity  and  the  pressure 
must  ce  continuous.  hence,  w ? th  (2,o) , 
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At  the  free  surface,  z •=  h,  a pressure  condition  hn=;  to  he  satis- 
fied. However,  we  are  only  interested  here  in  internal  waves 
whose  amplitude  nearly  vanishes  at  the  free  surface.  Therefore, 
it  is  permissible  to  write  for  the  boundary  condition  at  the  free 
surface 

w - 0 at  z - n 


To  consider  a particularly  simple  case  it  will  oe  assumed 
that  in  the  upper  layer  with  variable  density  p is  independent 
of  z , 


<p  = cor?  S'/' 


Since 
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it  follows  that  to  a very  close  decree  of  approximation 
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uhus  the  assumption  of  constant  p implies  a linear  distribution 
of  , the  quantity  which  is  actually  determined  in  ocear.ocraphic 
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where 


Since  in  the  lower  layer,  1 , (p  is  assumed  to  be  zero,  while 
it  is  constant  ir.  the  upper  layer,  2,  and  since  furthermore 
vanis.nes  at  z = 0,  w 0 at  z = h,  one  finds  for  the  vertical  velo- 
city from  (2.3 a ) r. ha t 

cuf  = /),  siwJ T^cr/k. & /)t \?icr sk, 

Uf^  = s2n  Tfi'ft-h) 

The  simplified  expression  for  w-^  follows  because  the  depth  of'  the 
fluid  is  much  smaller  than  the  wave  length  for  the  waves  considered 
here . 

At  the  internal  Doundary,  z = h^,  the  vertical  velocity  and 
its  derivative  with  respect  to  z must  be  continuous.  Hence, 
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and 
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Further,  the  vertical  disclacemcnt  C is  related  to  w by 

as 


= ur 


Hence 
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and  L}ie  amplitude  of  the  vertical  oscillation  is  at  the  internal 


boundary  z = hh  , 


C = ~ \(rh, 

L(T 


If  this  constant  is  used  as  tne  remaining  integration  constant 
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The  relation  between  wave  length  and  period  is  given  by  (2.9) 
which  is  a transcendental  equation  with  an  infinite  number  of  roots. 
It  can  be  solved  in  individual  cases  by  approxima tion  methods. 

Of  practical  importance  is  the  case  that  the  lower  layer  is 
much  thicker  than  the  upper  one.  in  this  case  a goon  approximation 
can  be  obtained  in  the  following  way.  Since  the  right- hand  side 
of  (2.9)  is  quj^e  large  and  negative 

7i  4 = +£ 

where  n is  an  integer  and  £ a small  quantity.  iionco,  approximately , 
If  this  expression  is  substituted  in  (2.9)  it  follows  that 


Since 


and  since 
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if  »<?*■  it  follows  that  the  wave  length 
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1.1 


As  a numerical  example  let 

= 2.1C"^  cm'1 

hj_  - 2700  m,  hp  - L'.O  m. 

Latitude  17°  GT  - 2 jr/l2.l\  hours 

The  figures  correspond  approximately  to  those  at  Meteor  Station 
385  (see  below),  except  for  the  surface  layer  where  the  density 
increase  is  more  rapid  than  is  assumed  here. 

For  lack  of  any  better  information  about  the  distribution  of 
the  internal  tidal  waves  perpcndi cular  to  the  direction  of  propa- 
gation, it  may  be  assumed  that  (X.  - 0.  Then  for  various  values 
of  n the  following  wave  lengths  woulu  have  to  exist  in  the  case 
of  Internal  waves  with  the  period  of  the  lunar  semidiurnal  tide 

n 0 1 2 3 

wave  length  ol  20  12  8.6  km 

From  the  foregoing  discussion  it  is  evident  again  as  from 
the  discussion  of  the  current  vector  that  for  a complete  check 
of  the  theory  simultaneous  observations  from  more  than  one  snip 
are  required,  in  oraer  to  determine  the  wave  length  and  the  dis- 
tribution of  the  oscillation  in  lateral  direction. 

If  such  internal  oscillations  occur  they  must  show  themselves 
also  in  periodic  changes  of  conservative  properties.  Ifxd  denotes 
such  a conservative  property, 


j.  n 

d"t  c>x 
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ferce  local  variations, 


, of  the  property  can  in  general  be 


due  either  to  horizontal  or  to  vertical  motions.  If  the  horizon- 
tal gradient  of  the  property  is  small  the  terms  depending  on  hori- 
zontal motion  can  be  neglected.  Further,  it  can  be  assumed  that 


where  the  first  term  on  the  right-hand  side  stands  for  the  mean 
value  of  at  each  level  z on  which  the  small  variation  " due 
to  the  internal  wave  is  superimposed.  Then 

— '—nr 


or  if 


HzSL  - — (jf  Q^£> 

crt  <3^- 

denotes  the  vertical  displacement  of  the  particle 


• § 


can  be  determined  from  the  variation  of  ^ and  the  ver- 


tical gradient  of  the  mean  value  o . 

In  the  lollowing  analyses  we  shall  confine  ourselves  to  the 
temperature  variations,  as  explained  previously.  in  the  case  of 
the  Atlantis  data  time  cross  sections  of  the  isotherms  were  already 
available  so  that  it  was  possible  to  analyze  the  heights  of  given 
isotherms  rather  than  the  temperatures  at  given  heights.  There- 
fore, the  same  method  has  been  used  in  the  analysis  of  the  Meteor 
data,  although  temperatures  at  given  heights  were  studied  in  the 
case  of  the  AlLair  data  in  order  to  save  time  and  labor.  It  need 
hardly  be  mentioned  tnat  the  temperature  is  not  a strictly  conser- 


n 


vativo  property  of  sea  water;  but  for  the  purposes  of  the  present 
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investigation  the  errors  arising  cut  of  this  assumption  are  cer- 
tainly not  greater  than  those  due  to  other  causes. 


ill.  Statistical  Procedure 

The  numerical  procedure  for  studying  the  various  series  of 
observations  ana  the  occurrence  or  non-occurrence  of  internal  oscil- 
lations has  been  described  briefly  in  the  preceding  report 
(Haurwitz,  1952).  Since  additional  comments  on  the  method  arc  to 
be  made  here  it  will  be  best  to  review  the  whole  procedure. 

If  a period  of  a certain  length  is  suspected  in  a set  of  data 
it  can  be  subjected  to  a harmonic  analysis  according  to  well-known 
methods.  As  a result  one  will  obtain  an  expression,  or  a sum  of 
expressions,  of  the  form 

sfft'  (3  #1) 


Here  t is  the  independent  variable,  for  instance  the  time,  expressed 
in  angular  measure,  n an  integer  representing  the  frequency,  and 
an  and  bn  the  harmonic  coefficients.  An  alternate  form  of  the 
foregoing  expression  is 


where 
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Since  it  is  always  possible  to  represent  a set  of  data  by  a 
series  of  terms  of  the  form  (3.1)  or  Its  alternate  (3.2)  the  har- 
monic analysis  does  not  prove  anything  about  the  reality  of  the 
suspected  periods.  r.ven  a set  of  random  numbers  can  be  represented 
by  a harmonic  series. 

In  many  cases  the  existence  of  certain  periods  in  the  observa- 
tion material  will,  of  course,  oc  highly  likely  because  of  the 
existence  of  a generating  force.  In  the  present  case,  for  instance, 
it  is  known  that  the  ocean  is  subjected  to  the  tidal  action  of  the 
sun  ana  tne  moon.  Therefore,  it  can  be  expected  that  the  oceano- 
graphic parameters  show  always  variations  with  the  tidal  periods, 
but  these  variations  may  be  so  minute  that  they  are  completely 
overshadowed  by  random  variations,  and  that  the  results  of  a har- 
monic analysis  represent  merely  the  effects  of  such  random  varia- 
tions. Hence  a statistical  test  is  necessary  to  decide  whether 
the  computed  harmonic  coefficients  represent  a real,  period  or  only 
random  fluctuations. 

Such  a test  for  reality  is  provided  by  the  expectancy  method. 
The  expectance  is  a measure  for  the  amplitude  An  of  a harmonic 
component  to  be  expected  in  ransom  data.  The  theory  of  probability 
shows  wbat  the  probability  is  that  the  actual  amplitude  exceeds 
the  expectance  by  a certain  amount.  Evidently  if  this  probability 
i?  very  small  it  is  unlikely  that  the  computed  amplitude  is  due  to 


random  data,  or  in  other  words, 
in  the  analyzed  data. 


that  the  period  exists 
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The  expectance  is  given  by  the  following  expression 

= 2/t^/V  n"1 

wher  is  the  standard  deviation  of  the  individual  n r.umbors 

(for  instance,  Chapman-Martels , 1940  or  Stumpff,  1940).  h-::-  would 

represent  the  mean  amplitude  in  such  a set  of  random  data.  it 
can  be  shewn  that  the  probability  that  any  amplitude  is  k times 
the  value  of  E-::-  is  given  by 


The  application  of  the  expectance  criterion  is,  however,  restricted 
to  those  cases  where  the  individual  members  of  the  set  of  numbers 
are  independent  from  each  other.  This  is  true  for  random  numbers, 
but  not  for  the  majority  of  geophysical  data.  In  the  cases  to  be 
considered  here,  for  instance,  the  subsequent  heights  of  the  iso- 
thermal surfaces  will  have  a curtain  autocorrelation  since  evi- 
dently large  positive  deviations  from  the  mean  value  are  more  likely 
to  be  followed  by  other  large  positive  deviations  from  the  mean  than 
by  large  negative  ones. 

It  is  possible  to  avoid  or  at  least  to  reduce  greatly  the 
effects  of  the  autocorrelation  by  introducing  instead  of  the  expec- 
tance a somewhat  different  expression,  the  expectancy.  In  order  to 
determine  this  quantity  let  the  whole  analysis  interval  tn  - t ^ 
be  an  integral  multiple  of  the  period  T whose  reality  is  to  be 
determined , 


h - *0  =S'  T 


Id  - 


1 


This  can  always  be  arranges  by  suitable  choice  of  the  end  points 
tp,  tQ.  The  total  interval  is  now  subdivided  into  equal  parts 
each  of  which  is  either  eoual  to  the  period  T or  to  a multiple 


of  T.  The  harmonic  coefficients  an 


and  fa- 


are  ohen  aeter- 


mined  separately  for  each  partial  interval  of  which  there  may  be 


r,  say,  wnere 


r=  . Let 

a(d2  = a<n2  , bd)2 


the  amplitude  for  each  partial  interval.  The  determination  of  the 
different  is  based  on  a number  of  points  which  with  the 

exception  of  the  last  points  of  one  partial  interval  and  the  first 
points  of  the  next  partial  interval  are  pretty  far  apart  so  that 
the  effects  of  any  autocorrelation  in  the  data  is  greatly  reduced. 
Therefore,  it  can  be  assumed  that  the  AnVJ-'  • determined  for  the 
different  partial  intervals  would  in  the  absence  of  the  suspected 
periodicity  behave  like  random  numbers  and  scatter  around  a mean- 
value  which  may  be  assumed  zero  a oriori , The  harmonic  coefficients 
ar/ ^ and  bn^^  can  be  consicerea  as  the  components  of  a vector 

which  may  be  plotted  in  a coordinate  system  with  ar  and  b_  as  the 
horizontal  and  vertical  axes,  respectively.  This  representation 
is  known  as  the  harmonic  dial,  wticn  win  be  discussed  later, 
when  the  various  a n v 1 ^ , b-,  ' - ' are  plotted  on  such  a r a rmoni o 

dial  the  end  points  of  the  vectors  will  form  a point  cloud  which 
should  in  the  case  of  complete  randomness  of  the  data  be  distributed 
around  the  origin  with  circular  symmetry. 

The  quadratic  mean  value  of  the  lengths  of  the  vectors  which 


represent  the  harmonic  coefficients  of  the  individual  partial 
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interve] s 


is  called  the  expectancy  of 

e2  = 2Z  A J‘)2/r 
i = 1 


the  se 


indi'vidual  vectors, 


Stumpff  1 V 4 0 ) has  shown  that  it  is  convenient  to  introduce  the 
expectancy  E {'or  the  distance  of  the  center  of  gravity  of  r.he 
cloud  of  r points  an^)  , br  ^ from  origin, 


The  distance  of  the  center  of  gravity  of  the  point  clouo  from  the 
origin  is  given  by 


(3.6) 


and  the  probability  that 


A„  = k 


(3.7) 


1 s 


(3.6; 


for  random  data.  Hence,  if  the  probability  p is  small  it  is 
likely  that  the  data  are  not  ranuomiy  distributed,  but  that  the 
suspected  period  is  real. 
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Since  the  derivations  which  lead  to  the  foregoing  formulae 
are  based  on  statistical  theory  they  presuppose  that  the  number 
of  intervals,  r,  on  which  the  test  is  based  is  not  too  small. 

This  can  be  seen  directly  if  the  following  case  is  considered. 

Suppose  that  the  an'^  and  the  for  all  partial 

intervals  are  the  same  and  that  the  amplitude  for  each  partial 
interval  is  An.  In  this  case  p = e“r.  If  r = 2,  for  Instance, 
the  probability  is  one  in  7 that  any  amplitudes,  even  though  they 
are  the  same  for  both  partial  intervals,  are  due  to  random  data. 
The  expectancy  test  is,  of  course,  not  applicable  to  so  few  data, 
and  the  example  merely  illustrates  this  point.  We  shall,  however, 
have  to  apply  it  later  to  data  where  r - 1;  and  where  therefore 
the  lowest  probability  is  one  in  that  the  amplitudes  found  by 
the  harmonic  analysis  are  aue  to  random  aata.  hven  in  this  c«.se 
the  application  of  the  test  can  only  be  justified  because  longer 
series  are  not  available. 

Another  point  to  be  discussed  here  concerns  the  exact  length 
of  the  perloa  wnose  realicy  is  to  be  tested.  If  a sufficiently 
long  series  of  observations  is  available  the  choice  of  a period 
slightly  aiflerent  from  the  exact  one  would  manifest  itself  in  a 
systematic  change  of  the  phase  constant  as  each  consecutive  period 
interval  is  analyzed.  hut  in  a short  series  comprising  only  a few 
periods  this  variation  of  the  phase  constant  would  be  completely 
obscured  by  the  irregular  changes  aue  to  superimposea  random  lluc- 
tuations  in  the  data. 
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In  thj  particular  problem  studied  here  this  difficulty  may 
arise  because  cf  the  slight  difference  between  the  period 
T0  = 12.42  hours  of  the  principal  lunar  tide  ana  the  period 
T = 12  hours  of  the  principal  solar  tide.  while  the  former  is 
as  a rule  the  dominant  one  this  is  not  always  true,  and  it  can 
therefore  not  definitely  be  ascertained  from  a short  series  of 
observations  which  of  the  two  periods  is  the  correct  one,  even 
if  the  expectancy  test  is  positive.  In  order  to  see  how  the 
choice  of  a period  slightly  different  from  the  correct  one,  say 
Te,  affects  the  oxpectaney  test  assume  that  a time  series  is  given 
by  the  simple  expression 

f(t)  = A cosy  t 

where 

V - 12  hr s/Te 


and  t is  the  mean  solar  time  in  angular  measure.  Let  the  partial- 
analysis  interval  be  T = 12  hours.  Then 


v = 1 - °-°34  = 1 -£ 

The  harmonic  coefficients  for  the  period  T computed  from  the 
(k  + 1)  th  anfll vsi « interval  are 


2ir(Jn-i) 

zrrJt 


TTXA\i) 

{' 

ETTA 


tcft  ■ 


z^ars;  z(/—£,)g> 

ire.  z-e 


&E22HS  JL . cl  s*n[Fg(zA+\)) 

ITS  2-£ 
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because  £ is  small  these  expressions  may  be  written  approximately 
as  follows 

oJt,  (/-§)a.CM [TZiUtij] 
J^~(/*§r)a  sulDts&Ah)] 

These  two  expressions  show  that  for  the  fictitious  period  T the 
amplitude  as  determined  from  each  partial  analysis  interval  remains 
nearly  constant  and  nearly  equal  to  A while  the  phase  constant 
changes  by  21(S  from  one  interval  to  the  next. 

If  the  coefficients  a^^  and  b^K  are  computed  for  the  first 
four  intervals  of  length  T one  finds  that 

p = 0.022 

while  the  lowest  value  for  the  probability  of  randomness  for  r = J4 
is  O.Olfl,  Evidently  . n i f'^'erer.c-.  between  thane  t*vu  values  is 

too  small  to  permit  any  decision  about  the  correctness  of  the 
assumed  period  length. 

but  even  with  a longer  series  of  observations  the  expectancy 
test  cannot  bn  used  for  a decision  about  the  correct  period  length. 
If  in  the  foregoing  example  11  consecutive  periods  of  the  length 
T were  used  one  would  find  that  p = 0.0011  while  the  theoretical 
minimum  value  is  0,000017.  The  difference  is  here  considerably 
larger.  However,  since  in  actual  cases  the  observations  would 
also  be  subjected  to  random  effects  a value  of  p larger  than  the 
theoretical  minimum  can  hardly  be.  attributed  to  the  wrong  choice 
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of  period.  On  the  other  hand,  with  such  e.  longer  series  the  indi- 
vidual harmonic  coefficients  wou.lu  show,  despite  superimposed  ran- 
dom effects,  a systematic  change  of  the  phase  corstant  if  the 
assumed  period  deviates  r rom  the  actual  one  present  in  the  obser- 
vation material. 

'The  foregoing  remarks  require  some  modifications  if  the  par- 
tial analysis  intervals  consist  of  more  than  one  period.  But  since 
the  series  analyzed  here  are  not  long  enough  to  permit  such  a com- 
bination of  the  data  this  case  need  not  be  considered  here. 


IV.  Analysis  of  the  Data 
A summary  of  the  stations  which  were 
Table  I,  Station  3^91  of  Atlantis  Cruise 

Table  I 

LIST  Of  STATIONS 


Investign  ted 
which  was 


77 


Designation 
Atlantis  Cruise  77 
Meteor  Station  385 
Meteor  Station  J;38 
Altair 


Location 

34 « 1°N  o 5 . 9 ° w 
16. 6° N 46. 3° 

30.0°N  U3-8°w 
LJi.5°N  lit . 0 0 a' 


Duration 
16-29  June  1938 
12-1U  Feb.  1936 
26-28  Apr.  1938 
16-20  .Tune  193-8 


is  given  in 
discussed 


Reference 
Kaurwitz,  1952 
Schubert,  19kb 
Schubert,  1944 
Do  f ant.  t 9q.Ca 


In  the  earlier  report  has  been  included  for  completeness.  The 
average  vertical  density  distributions  as  obtained  from  observa- 
tions made  at  these  stations  are  shown  in  Fig.  1.  In  the  case  of 
the  Atlantis  and  the  Meteor  observations  one  series  of  measure- 
ments was  extended  to  much  greater*  depths.  The  results  of  these 


deep  soundings  are  indicated  by  individual  points  not  connected  by 
a continuous  curve.  In  the  upper  layers  where  mean  values  frcm 
the  serial  observations  are  also  available  the  observations  of 
the  single  deep  sounding  agree  very  woh  with  the  means  of  the 
serial  observations.  The  average  depth  during  each  station  is 
indicated  by  the  horizontal  line  at  the  bottom  of  eacn  diagram. 

The  number  gives  the  depth  in  meters.  in  the  case  of  the  Altair 
which  anchored  over  a seamount  the  greatest  and  smallest  depths 
are  shown. 

The  vertical  density  distribution  during  the  Atlantis  station 
shows  a very  rapid  density  increase  downward  in  the  uppermost 
100  meters,  followed  by  a more  homogeneous  layer.  Between  50C 
and  1200  m the  density  gradient  is  again  slightly  larger.  Simi- 
larly, the  Meteor  Station  385  and  the  Altair  Station  show  a surface 
layer  of  large  density  gradient  followed  by  a more  homogeneous 
layer.  At  Meteor  Station  385  the  density  at  the  sea  surface  is 

r\  j r n "1  t r*  V>  A 1 »»  V-.  i — • l—  »>  — ■>  ■*- 1 — - •»  ^ « T T a • i V-  ■ • — — — • — • 
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below  it.  \t  Meteor  Station  ip 38  the  density  gradient  in  the  highest 
layers  is  much  less  strong  than  at  the  other  stations  and  con- 
siderably mnw  uniform  through  the  first  1000  m.  The  temperature 
distribution  resembles  the  distribution  of  the  density  very  closely 
In  all  cases. 

e.  Review  of  Atlantis  Data 

In  order  to  recapitulate  briefly  the  results  of  the  stuuy  of 
the  Atlantis  data  and  to  illustrate  the  method  of  presentation  of 
the  results  of  the  present  investigation,  Fig.  2 shows  the  harmonic 
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a. 


dial  for  the  15°C  Isotherm.  The  average  depth  of  this  isotherm 
was  630  m.  The  analysis  of  the  total  time  interval  of  observa- 
tions which  composed  21  lunar  semiaiurnal  tidal  periods  was  made 
for  seven  groups,  each  consisting  of  3 consecutive  tidal  periods. 
The  points  denoted  by  A,  B,  . G represent  the  results  of  the 
individual  analyses.  The  mean  of  all  21  periods  is  represented 
by  the  cross  connected  by  a straight  line  to  the  origin  of  the 
coordinate  system^.  This  mean  amplitude  A.,  is  only  1.21  times 
larger  than  its  expectancy.  Hence  the  probability  that  it  is  due 
to  random  data  is  0.22,  far  too  large  to  regard  the  computed  ampli- 
tude for  the  lunar  tidal  period  in  the  data  as  real.  Instead  of 
showing  the  expectancy  of  Am  in  P'ig.  2 circles  of  equal  probability 
have  been  drawn  with  the  origin  as  the  center.  The  radius  of  each 
of  these  circles  indicates  the  probability  that  a value  of  Am  of 
the  same  magnitude  is  due  to  random  data. 

The  results  of  the  analyses  of  the  individual  7 groups,  repre- 
sented h v th°  points  A through  C,  scatter  widely.  In  oi'der  so  nave 
a measure  of  this  scatter  the  deviations  ^ a^^  and  /S  ^ of  the 


individual  ali;  and  b. 


'rcir.  the  nean  have  been  formed.  Then  the 


quantity 


has  been  computed.  The  resulting  number  nay  be  considered  as  the 


* In  Table  I of  the  previous  report  the  mean  value  of  the  bg  for 
the  15°  isotherm  should  be  -l.o3  m,  rather  than  -1.9C  m.  Conse- 
quently in  Table  II  of  that  report  Am  = k = 1.23,  and  p = 

0.22,  a change  of  no  significance. 
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probable  error  uf  A^.  Its  magnitude,  5*2  in,  is  indicated  by  the 
circle  around  the  end  point  of  as  the  center.  The  fact  that 
the  probable  error  is  larger  than  the  value  of  A,,,  shows  again 
that  no  tidal  period  is  present  in  the  Atlantis  data. 

The  harmonic  dial  shows  also,  on  the  margin  of  the  figure, 
for  each  position  of  the  harmonic  coefficient  that  time,  in  loca?. 
lunar  hours,  at  which,  the  maximum'  would  occur.  because  of  our 
choice  of  the  expression  ror  the  harmonic  representation,  (3.2), 
the  hours  Increase  in  the  counterclockwise  direction,  opposite  to 
the  conventional  clock  faces. 

b.  Meteor  Station  385 

The  location  of  this  station  is  given  in  'iable  I.  The  serial 
observations  extend  only  through  the  uppermost  150  m.  Therefore 
the  analysis  has  to  be  confined  to  thi^  layer.  As  already  pointed 
out  by  von  Schubert  (1944)  an  Inspection  of  the  figures  showing  the 
heights  of  isothermal  and  isohaline  surfaces  as  functions  of  time 
(his  Figures  7 and  8)  shows  fairly  clearly  a semidiurnal  tidal 
Deriod.  Hence  it  could  be  surmised  that  the  expectancy  teso  should 
in  the  present  case  give  a positive  result  as  to  the  reality  of  the 
semidiurnal  lunar  tide.  In  orcer  r.n  ncr>fnr*m  nhp  nnnl«ci ; time 

A \J 

cross  sections  of  the  Isothermal  surfaces  24.5°^,  23°C,  2?°C,  21°C. 
and  20°C  were  drawn  on  the  basis  of  the  published  soundings,  and 
values  of  the  heights  of  these  surfaces  were  read  off  for  each 
lunar  hour  beginning  with  C100  local  lunar  time,  corresponding  to 
12  February  1938,  11:15  45°w  meridian  time.  Altogether,  hourly 
sounuings  were  taken  lor  60  solar  hours,  but  because  ol  the  greater 
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length  of  the  lunar  hour  only  four  semidiurnal  lunar  tidal  periods 
were  available.  The  harmonic  analysis  was  performed  separately 
for  each  individual  tidal  period  and  the  results  are  listed  in 
Taole  II.  The  table  also  gives,  in  the  lowest  line,  the  average 
depth  of  the  various  isothermal  surfaces  during  the  observation 
period  , 

From  these  harmonic  coefficients  the  mean  amplitude  A^,  its 
expectancy  E,  the  phase  angle the  probability  p,  that  Am  is 
oue  to  random  data  and  the  probable  error  of  Am,  P.E.  (Am)  have 
been  computed.  They  are  given  in  Table  III.  The  main  conclusion 
to  be  drawn  from  this  table  is  that  it  is  not  very  probable  that 
the  amplitudes  of  the  semidiurnal  oscillation  as  computed  are  due 
to  random  data,  in  other  words  thai.  the  exlstance  of  semidiurnal 
internal  tide  is  very  probable.  The  value  of  the  probability  p 
that  the  amplitude  Am  is  due  to  random  data  is  between  0.02  and 
0,03  for  the  different  Isothermal  surfaces.  This  probability  for 
randomness  is  not  very  small,  but  it  should  be  remembered  that  with 
only  four  groups  of  data  the  lowest  probability  which  the  expectancy 
test  could  give  is  e-^  = 0,018. 

Table  III 

EXPECTANCY  TEST  FOR  THE  LUNAR  SEMIDIURNAL  TIDE 


METEOR  STATION  365 


ISOTHERM 

23°  C 

22°  C 

21°  C 

ru 

0 

0 

0 

Am 

1 1 . c.  m 

13 . 5 

14*4 

13.9 

1 1 r* 

- J-  • ^ 

242° 

249° 

25  0° 

251° 

257° 

E 

5.7  m 

o.9 

7 hi 

7.5 

6.3 

P 

0.02 

0.02 

C.G2 

0.03 

0.03 

F.E, ( A-) 

1.5  m 

: j 

}- 1 

— 

HI 

j 

2 * * 

V, 


1 
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For  the  2jcC  isotherm  a harmonic  dial  similar  to  Fig.  2 for 
the  Atlantis  data  is  shown  in  Fig.  3.  A comparison  of  the  two 
figures  shows  clearly  hew  much  more  probable  l he  existance  of  the 
semidiurnal  lunar  tide  is  in  the  data  from  the  Meteor  Station  3$5 
than  in  the  Atlantis  data. 

At  smaller  depths  than  about  80  m,  the  mean  depth  of  the  24.53C 
isotherm,  tidal  oscillations  could  not  be  observed,  as  has  been 
pointed  out  already  by  von  Schubert  (1944)-  Sven  if  tidal  oscil- 
lations were  in  existence  here  it  would  be  very  difficult  to  observe 
them  because  the  vertical  gradients  of  temperature  and  salinity  are 
very  small. 

As  Table  III  shews  the  amplitude  increases  downwards  and  reaches 
a maximum  tor  the  22°C  isotherm  while  the  phase  angle  Increases 
steadily,  although  somewhat  irregularly.  But  the  differences  are 
not  larger’  than  the  probable  errors.  hence,  it  must  be  concluded 
that  no  significant  variations  with  depth  of  the  internal  tidal 
oscillation  can  be  found.  Since  the  layer  from  which  data  are 
available  comprises  only  $0  ir.  this  is  not  surprising. 

As  pointed  out  in  section  3 It  1-  not  possible  to  decide  on 
the  basis  of  the  3hort  series  of  available  data  whether  the  period 
length  is  really  equal  to  that  of  the  principal  lunar  tide.  In 
order  to  illustrate  this  difficulty  on  actual  data  a harmonic  ana- 
lysis was  performed  with  the  principal  solar  tidal  period  as  the 
analysis  interval.  For  this  Analysis  not  the  height  of  an  isothermal 
surface  was  chosen,  but  the  temperature  variation  at  120  m depth 
which  could  be  taken  directly  from  von  Schubert's  (1944)  tables. 
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Since  the  whole  series  comprises  6C  hours  five  individual  periods 
are  available  for  the  expectancy  test.  The  probability  that  the 
mean  amplitude  at  this  level,  C.97°C  is  due  tc  random  fluctuations 
is  0.02.  This  figure  is  slightly  smaller  than  the  probability, 

0.03  that  the  lunar-tidal  harmonic  coefficients  of  the  21°C  end 
20°C  isotherms  are  due  to  random  f 1 nc tuat i -ns . But  because  of  the 
short  length  of  the  series  of  observations  the  slight  difference 
between  these  probabilities  has  no  significance. 

In  order  to  demonstrate  that,  on  the  other  hand,  the  expectancy 
test  shows  the  nonexistence  of  periods  whose  lengths  differ  sub- 
stantially from  the  real  ones  the  temperature  variations  at  120  m 
depths  were  also  analyzed  for  a period  of  eight  hours.  The  expec- 
tancy test  showed  that  the  resulting  mean  amplitude  of  0.1):  °C  would 
occur  in  random  data  with  the  probability  0.79.  Hence  such  a period 
does  not  exist  in  the  data. 

c.  Meteor  Station  438 

The  exact  location  of  this  station  which  is  about  13°  north  of 
Meteor  Station  385  is  shown  in  Table  1.  As  Fig.  1 shows  there  is 
no  layer  of  very  rapid  density  increase  near  the  surface.  hut  a 
pronounced  density  gradient  e?,-.  t t-fs  througnout  the  highest  800  m so 
which  depth  the  regular  serial  observations  were  made.  As  at  the 
other  Meteor  Station  observations  were  made  at  hourly  intervals, 
from  26  April  193<3,  10:12  45°*  meridian  time  to  28  April  1938, 

20:58  45 meridian  time,  a period  of  almost  60  hours.  Thus  again 
four  consecutive  lunar  semidiurnal  periods  could  be  analyzed.  As 
?ero  for  t 1400  local  lunar  time  or:  26  April  was  chosen  which 
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corresponds  to  10:4b  450|V  meridian  time.  The  temperature  data  pub- 
lished by  von  Schubert  (1944)  i-iere  again  used  to  construct  a time 
cross  section  of  the  height  of  various  isothermal  surfaces,  viz., 
18.5°,  l8°C,  17.5°C,  17 °C , etc.,  to  10°C  from  which  the  height  of 
each  surface  for  every  lunar  hour  could  be  reaa.  -ven  a casual 
inspection  of  this  time  cross  section  showed  that  internal  oscil- 
lations of  tidal  character  if  present  at  all  are  much  less  well 
developed  than  at  Meteor  Station  385a 

Table  IV  shows  the  harmonic  coefficients  for  the  four  partial 
intervals.  The  results  of  the  expectancy  test  are  summarized  in 
Tcble  V whose  arrangement  is  similar  to  Table  ill,  but  probable 
errors  for  the  Am  and  the  mean  phase  angles  have  not  been  computed. 
Except  for  the  16°C  and  14°C  isotherms  the  probability  is  always 
greater  than  one  half  that  the  amplitudes  found  by  harmonic  analysis 
are  due  to  random  fluctuations.  Even  for  the  16°C  and  14°C  iso- 
therms the  probability  of  randomness  is  not  small  enough  to  regard 
the  existence  of  a semidiurnal  oscillation  as  established. 

Table  V 

EXPECTANCY  TEST  FOR  THE  LUNAR  SEMIDiuHNAL  TIDE 


i-u-i  j.  li'jn  oirtiiww  4,50 


ISOTHERM 

18, 9°C 

1 «-7  O r* 

1 I • *2  ^ 

16°C 

ll|°C 

H 

l\J 

o 

O 

o 

o 

Am 

3.1m 

1.0 

5.0 

1.9 

1.9 

2.4 

E 

4 . 0 m 

| 

1.5 

3.3 

1.4 

4.9 

3.3 

P 

°-55  i 

0.63 

0.11 

| C » 16 

0.65 

0,58 
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It  may  t3  pointed  out  that  the  mean  vertical  temperature  gra- 
dient has  the  following  values: 

I 

Depth  0 - 100  n luO  - 300  m 300  - qOO  ra  400  - 800  n: 

°C/10Cm  .91  04  1.03  1.58 

Hence  the  local  temperature  changes  from  which  the  height  changes 
of  the  isothermal  surfaces  have  been  computed  arc  very  small;  even 
the  amplitude  of  the  lb°C  isotherm,  5 m,  corresponds  only  to  about 
0.07°C.  The  appearance  of  greater  randomness  of  the  higher  iso- 
therms, lo.5cC  and  1 7 . 5C C , might  therefore  be  attributed  to  th6 
still  smaller  lapse  rates  in  these  layers.  However,  if  this  reason- 
ing were  correct  the  1?°C  and  10°C  isotherms  should  also  show  a 
lower  probability  of  randomness  than  they  actually  ao  since  they 
are  in  the  layer  between  400  and  800  m where  the  temperature 
decreases  even  more  rapidly  than  between  300  and  400  m.  it  is,  of 
course,  possible  that  the  existence  cf  an  Internal  tidal  oscillation 
was  confined  to  the  vicinity  of  the  16°C  isotherm  because  of  its 
position  near  a discontinuity  of  the  temperature  lapse  rate.  but 
because  of  the  rather  large  value  of  the  probability  p that  the 
amplitudes  are  due  to  random  data  rather  than  an  internal  tide  it 
does  not  appear  necessary  to  consider  this  possibility  here  In 
greater  detail. 

d.  Altair  Station  (1938) 

The  third  series  of  observations  was  obtained  by  the  Altair 
during  the  period  from  lo  to  20  June  1938  (Defant,  1940a).  Hourly 
soundings  were  made  for  a period  of  o7  hours  sc  that  7 complete 
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lunar  semidiurnal  tidal  periods  were  available.  The  hourly  sound- 
ings were  made  to  a depth  of  600  m.  Instead  of  determining  and 
analyzing  the  depths  of  various  isothermal  sur faces  as  in  the 
f'ese  of  the  Meteor  data  the  temperature  changes  at  given  levels 
will  here  be  used  directly  follov.’ing  the  procedure  used  by  Defant 
(1940b).  Further,  Instead  of  analyzing  separately  the  temperature 
data  at  certain  levels  the  temperature  of  the  following  levels  have 
boon  added  and  their  sums  analyzed: 


I 

25  m, 

50  m, 

7 5 m 

II 

100  m, 

150  m, 

200  m 

III 

300  m, 

400  m. 

600  m, 

800  m 

This  combination  of  data  has  the  advantage  that  it  smooths  out  3ome 
of  the  irregularities  in  the  temperature  variations  at  the  indivi- 
dual levels.  As  pointed  out  by  Defant  the  indicated  selection  of 
heignt  combinations  follows  the  actual  stratification  of  the  ocean 
at  the  Altalr  station  into  a surface  layer  down  to  100  m with  a 
sharp  discontinuity  at  the  surface,  a middle  layer  of  normal  den- 
sity Increase  to  about  350  n and  a lower  layer  with  normal  density 
increase  separated  from  the  middle  layer  by  a weaker  discontinuity 
layer. 

For  the  analysis  of  the  Altair  data  according  to  the  lunar 
semidiurnal  tidal  period  lunar  time  for  Greenwich  rather  than  for 
the  local  meridian  was  used.  The  origin  of  the  time  a::is  was 
09 : 00  lunar  Greenwich  time  corresponding  to  17  June  1938,  00:02  GMT. 


Ihe  results  of  the  analysis  for  the  three  layers  are  summarized 
in  Table  VI,  The  individual  and  the  mean  temperature  amplitudes 

Table  VI 

HARMONIC  ANALYSIS  FOH  THE  LUNAR  SEMIDIURNAL  TIDE 

ALTAIR  STATION 

16-20  June  1938  4/;.5°N  3u.0°* 


LAYER  I II  III 


Group 

£2 

— 

b 2" 

a2 

b2"" 

a2*' 

b2  ' 

A 

0.34 

-0.35 

0.15 

-0.19 

c.34 

-0.20 

B 

-0.16 

-0.18 

0.16 

-0.06 

0.02 

-0.09 

U 

-0,22 

-0.06 

-0,02 

-0.10 

0.02 

-0.01 

D 

0.30 

-0.40 

0.20 

-0.02 

0.10 

0.14 

E 

-0.08 

-0.20 

0.08 

0.06 

-0.05 

0.01 

F 

o.35 

0.16 

0.06 

0.18 

0.14 

-0.11 

G 

0.28 

0 . 64 

0.10 

-0.11 

0.23 

-0.24 

Mean 

0.  lib 

-0.056 

0.104 

-0.034 

0 , 114 

-0.057 

Mean  Temp, 

14. 

85°c 

12 

. 79°C 

IOji 

Ul£ 

The  units  of  the  harmonic  coefficients  for  layers  I and  II 
are  0.33°C,  for  layer  111  0.25cC. 

are  quite  small*  only  of  the  order  ui  a i'ew  hundredth  degree  Centi- 
grade. In  order  to  show  the  corresponding  magnitude  of  the  vertical 
displacements  the  mean  vertical  temperature  gradients  must  be  con- 
sidered. These  are  as  follows  according  to  Defant  (1940d): 

Depth  m 0 25  50  75  100  200  300  400  600  8CC 

°C/lC0m  3.9  5.3  4 . 2 1.6  1,0  1.2  l.o  1.5  0.7 


For 
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instance  since  the  mean  amplitude  for  layer  11  is  about  Oji_<o  = u 
the  corresponding  amplitude  of  the  vertical  oscillation  would  be 
less  than  ii  meters.  This  value  is,  of  course,  computed  on  the 
assumption  Lha t the  temperature  variation  h entirely  due  to  ver- 
tical motion  anu  net  to  horizontal  advecticn.  The  latter  possi- 
bility has  been  considered  by  Defant. 

The  expectancy  test  is  summarized  in  Table  VII. 

Table  VII 

EXPECTANCY  TEST  FOR  THE  LUNAR  SEMIDIURNAL  TIDE 
ALTAIR  STATION,  3958 


LAYER  I 

ii 

III 

Am  0.043°C 

0.036 

0.032 

k i o.o54°c 

0.022 

0,021 

p 1 0.53  (0.58) 

0.06  (0.03) 

. . _ 

0.10  (0.12) 


The  probability  is  very  high  that  the  amplitude  for  the  semi- 
diurnal lunar  tide  computed  for  the  top  layer  is  not  real.  The 
probabilities  for  the  middle  and  lower  layers  are  considerably 
smaller  giving  greater  confidence  in  the  reality  of  the  Internal 
tide  at  these  levels,  even  thougn  tr.e  values  of  p are  not  small 
enough  to  make  the  existence  of  the  internal  tide  here  a certainty. 
But  despite  the  relatively  high  values  of  p for  these  two  layers 
it  l=i  remarkable  that  especially  in  layer  II  where  the  temperature 
gradients  are  smallest  the  probability  of  randomness  is  also 
smallest.  This  fact  suggests  strongly  that  with  a longer  series 
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* 

V 

1 


L 
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z 
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t 
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\ 

r" 
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of  observations  tne  existence  of  the  internal  semidiurnal  tide 
would  have  been  brought  out  by  the  expectancy  test.  The  numbers, 
p,  in  parentheses  will  be  discussed,  later,  un  page  39.  1 

From  an  inspection  of  the  current  measurments,  and  the  tem- 
perature and  salinity  soundings  Defant  (19i-'f0b)  surmised  that  a 
17-houriy  wave  was  superimposed  c:i  the  semidiurnal  lunar  tidal 
oscillation.  That  an  oscillation  of  1?  hours  duration  might  exist 
is  theoretically  plausible  slne^  the  period  of  the  inertia  motion 
at  this  latitude  is  1?  solar  hours.  Therefore  the  Altair  oDserva- 
tions  have  also  been  analyzed  with  respect  to  this  period.  Five 
complete  inertia  periods  are  contained  in  the  time  Interval  during 
which  soundings  were  made.  For  the  origin  of  the  time  axis,  the 
same  time  was  chosen  as  in  the  case  of  the  analysis  for  the  lunar 
tidal  period}  viz.,  1 ? June,  00:02  GMT.  The  harmonic  coefficients 
are  given  in  Table  VIII 


Table  VIII 

HARMONIC  ANALYSIS  FOR  THE  IN ART! A PERIOD  {17  Hours) 

ALTAIR  STATION,  1938 

LAYER  I II ill 


Group 

a" 

b::’ 

a'" 

b';:' 

a"' 

b“ 

A 

0.11 

-0.36 

0.GU 

! -0.21 

0.22 

-O.lJi 

B 

0.11 

0.09 

-0.05 

-0.03 

-C.01 

0.09* 

C 

1 

0.27 

1 

-0.17 

0.10 

C.09 

-0.14 

C . 20 

D 

-0.30 

-0.17 

0.13 

-o.oB 

-0.C7 

-0.0b 

£ 

-0.R5 

0.12 

-0.02 

-0.03 

0.01 

0 . 1L 

Mean 

-0.053 


-0.10 

0.037 

i 

_ n r.£  i 1 

1 

0.000 

0.046 

The  units  of  the  harmonic  coefficients  for  layers  I and  II 
arc  C.33°C,  for  layer  III  0 .25cC. 
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The  expectancy  test  is  summarized  in  Table  IX.  Since  the  pro- 
bability is  large  that  amplitudes  of  the  magnitudes  found  for  the 
inertia  period  would  be  found  in  random  data  the  reality  of  the 
17  hour  period  is  open  to  serious  doubt. 

Table  IX 

EXPECTANCY  TEST  FOR  INERTIA  PERIOD 
ALTA1R  STATION,  1938 


LAYER 

I 

11 

III 

^-m 

0.03&°C 

0.021 

0.011 

E 

o.oi+7°c 

0.020 

0.019 

P 

0.53 

0.3U 

0.70 

In  order  to  test  how  well  the  superposition  of  the  two  periods, 
namely,  the  lunar  semidiurnal  and  the  inertia  period,  represents 
the  observed  data  tne  mean  temperature  of  layer  il  has  been  com- 
puted by  harmonic  synthesis  for  every  lunar  hour  and  compared  with 
the  observed  mean  temperature.  The  correlation  coefficient  between 
the  83  computed  and  observed  values  is  only  0.5U  indicating  that 
the  representation  is  not  very  satisfactory. 

The  existence  of  a 17-hour  period  in  the  Altair  data  appears 
thus  unlikely.  Nevertheless  the  possibility  of  the  existence  of 
two  superimposed  periods  in  the  observation  material  raises  the 
question  as  to  the  influence  of  the  coexistence  of  two  or  more 
periods  on  the  expectancy  test. 


V. The  Superposition  of  Two  Oscillations 


Suppose  that  a function  of  time  p(t)  is  the  sum  of  two  oscil- 


lotions  of  different  periods.  and  7,  , so  that 


rv  • \ 

^ i ) 


Jt  : b sin  -0(oo 


t <■  ^ sir 


r.  tr, 


where  £ - 2 77 i = £TT/  T^ . In  order  to  perform  the  expec- 
tancy test  with  respect  to  the  period  T one  would  determine  the 
following  coefficients  by  harmonic  analysis 


= -S 


'( k + 1 ) T 


-J  ^p(t)  cos  ^ t o t,  b = ^ j>(  t) 


J IcT 


(k+l)T 


rl^"t  a 


which  become,  on  performing  the  integration 


2 1 1 / r"  -1-/  -i  r — / 

A + ^^7tT2 j(Xcos[(  k+l/2 ) 2TT^-J+psin[(k+l/2)  2TT^-Ji  (i;.?) 

1 Kl  ^ L 1 1 J 


-rr- 

sinTT 


b = B + - 


sin  77“  ^ _ "*) 

2 Ti  J T T/7  t / 

+ —?VXsin[  (k  + l/2)  2JT  Scos[{k  + l/2)Zfr^  ] / 

1 " i^K  1 J 


v u . 3 ) 


If  the  expectancy  test  were  performed  not  by  dividing  the  whole 
analysis  interval  Into  groups  of  the  length  of  a single  period, 


T‘.«v* 


but  into  groups  which  comprise  an  integral  number  of  periods,  the 
foregoing  formulae  would  only  have  to  be  modified  slightly. 

The  last  formulae  for  a and  b show  that  because  of  the  pre- 
sence of  an  additional  oscillation  the  harmonic  analysis  will  in 
general  not  give  the  correct  values  A and  3 for  the  coefficients 
of  the  period  , but  different  values,  the  difference  depending 
on  the  ratio  of  the  periods  T ^ and  and  on  the  magnitude  of  the 
coefficients  0(  ana  ^ of  the  superimposed  period,  if  the  obser- 
vation series  is  long  enough  and  if  the  ratio  T /)  to  T,  is  an 
lnteiger  it  is,  of  course,  possible  to  eliminate  the  effect  of  a 
superimposed  period  when  performing  the  expectancy  test.  For 
instance,  if  = 24  hours  and  the  period  to  be  studied  = 12 

hours,  it  is  merely  necessary  to  subdivide  the  whole  set  of  data 
into  groups  of  the  length  2 * r"  x 12  hours,  in  order  to  eliminate 
the  effect  of  the  24-hourly  wave  on  the  expectancy  test. 

In  the  cases  which  have  been  studied  here  it  might  in  general 
be  suspected  that  a wave  of  approximately  2 4 hours  period  also 
exists  because  of  the  corresponding  terms  in  the  tidal  potentials 
of  the  moon  and  of  the  sun.  It  is  not  possible  to  eliminate  the 
possible  effect  of  such  a longer  period  here  by  the  combination 
of  a number  of  periods  since  the  total  lengths  of  the  series  are 
not  large  enough.  However,  if  an  appreciable  24-hou\ ly  oscilla- 
tion were  present  it  would  shotf  itself  in  systematic  variations 
of  the  coefficients  a£  ana  \>2  for  each  individual  time  interval 
of  12  lunar  hours,  bir.ee  no  such  variations  are  apparent  in  the 
computed  coefficients  it  must  be  concluded  that  the  expectancy 
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teste  have  not  been  influenced  by  an  oscillation  wiln  a period 
near  ?Ji  hour's. 

In  the  case  of  the  Altair  data  with  trie  suspected  inertia 
period  of  17  hours  matters  are  more  complicated  because  of  the 
non-integral  raeio  c:  the  two  periods  involved.  It  is  possible 
to  apply  a correction,  for  instance  to  the  harmonic  coefficients 
or  the  lunar  semidiurnal  tine  given  in  Table  VI,  on  the  basis  of 
equations  Up. 2)  and  (Ip.  3).  The  coefficients  a and  b may  be  identi- 
fied with  the  coefficients  as  listed  in  Table  VI.  The  correct  co- 
efficients would  be  A and  5 which  can  be  computed  if  the  coef- 
ficients <X  and  p of  the  17-hour  oscillation  arc  known.  To  make 
this  correction  these  coefficients  were  identified  with  the  mean 
values  of  the  coefficients  a and  d cf  the  17-hour  oscillation  listed 
* n Table  vru.  Such  a correction  has  been  applied.  The  corrected 
values  of  a2  and  b2  are  not  reproduced  here.  Only  the  changes  in 
the  probabilities  of  randomness  due  to  the  corrections  are  shown, 
namely  by  the  figures  in  parentheses  given  in  the  last  line  of 
Table  VII.  In  general  the  changes  due  to  the  corrections  are 
slight.  Only  in  the  case  cf  the  middle  layer  II,  the  probability 
that  the  computed  amplitudes  are  due  to  random  fluctuations  is 
halved  a fact,  which  lends  further  weight  to  our  previous  remarks 
about  the  reality  cf  internal  tide  in  this  layer.  It  need  hardly 
be  pointed  out  that  the  whole  correction  method  is  very  cruae  since 
the  time  variation  of  the  data  is  only  very  inadequately  repre- 
sented by  the  superposition  of  a lunar  semidiurnal  and  an  inertia 
oscillation.  The  correction  has  been  applied  here  mainly  to  antici- 
pate  a possible  objection  to  tne  expectancy  test. 
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VI.  Conclusion 

The  preceding  analyses  show  that,  internal  tidal  waves  appeared 
with  reasonable  certainty  only  in  one  case,  namely,  at  Meteor  Sta- 
tion 3q5,  while  at  the  other  stations  the  expectancy  test  gives 
either  inconclusive  results  or  shows  that  no  such  oscillations 
occurred. 


It  is  apparent  from  the  discussion  presented  here  that  only 
longer  series  of  data  lasting  more  than  just  2 or  3 days  permit  a 
decision  as  to  the  reality  of  internal  tidal  waves,  and  it  is  one 
of  the  principal  aims  of  this  report  to  stress  the  need  for  such 
adequate  observations . Furthermore,  the  value  of  any  series  of 
observations  will  be  greatly  enhanced  if  measurements  \re  taken 
simultaneously  at  three  observation  points,  as  explained  in  Section 
II.  Such  synoptic  observations  not  only  support  each  other, 

but  at  the  same  time  give  much  needed  information  about  the  wave 
length  and  direction  of  travel  of  these  internal  ««v>'-s  and  about 
their  lateral  extent. 

The  question  whether  these  internal  wsves  are  due  to  hori- 
zontal or  vertical  motions  or  to  a combination  of  both  has  only 
briefly  oeen  touched  here  since  it  ic  beyond  th*»  scope  of  t::is 


reoert . Also  for  s decision 
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three  points  would  be  very  helpful  since  they  would  give  the  neces 
sary  information  about  the  horizontal  gradients  of  the  required 
oceanographic  parameters.  In  addition,  the  observations  should 
also  include  cun  ent  observations  which  woulo  show  wnether  the 


horizontal  advection  together  with  the  horizontal  grauients  couia 


Moreover,  since 
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alone  account  for  the  observed  changes  or  not.  Moreover,  since 
the  internal  waves  must  produce  current  fluctuations  data  on 
these  currents  would  further  corroborate  any  tidal  fluctuations 
observed  in  temperature  of  salinity. 
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